Mon. Not. R. Astron. Soc. 000, \TM (2008) Printed 7 July 2008 (MN WT&L style file v2.2) 

Can a primordial magnetic field originate large-scale 
anomalies in WMAP data? 
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ABSTRACT 

Several accurate analyses of the CMB temperature maps from the Wilkinson Mi- 
crowave Anisotropy Probe (WMAP) have revealed a set of anomalous results, at large 
angular scales, that appears inconsistent with the statistical isotropy expected in the 
concordance cosmological model ACDM. Because these anomalies seem to indicate 
a preferred direction in the space, here we investigate the signatures that a primor- 
dial magnetic field, possibly present in the photon-baryon fluid during the decoupling 
era, could have produced in the large-angle modes of the observed CMB temperature 
fluctuations maps. To study these imprints we simulate Monte Carlo CMB maps, 
which are statistically anisotropic due to the correlations between CMB multipoles 
induced by the magnetic field. Our analyses reveal the presence of the North-South 
angular correlations asymmetry phenomenon in these Monte Carlo maps, and we use 
these information to establish the statistical significance of such phenomenon observed 
in WMAP maps. Moreover, because a magnetic field produces planarity in the low- 
order CMB multipoles, where the planes are perpendicular to the preferred direction 
defined by the magnetic field, we investigate the possibility that two CMB anomalous 
phenomena, namely the North-South asymmetry and the quadrupole-octopole planes 
alignment, could have a common origin. Our results, for large-angles, show that the 
correlations between low-order CMB multipoles introduced by a sufficiently intense 
magnetic field, can reproduce some of the large-angle anisotropic features mapped in 
WMAP data. We also reconfirm, at more than 95% CL, the existence of a North-South 
power asymmetry in the WMAP five-year data. 
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1 INTRODUCTION 

The five-year data from the 
Anisotropy Probe (WMAP ) 



Wilkinson Mi cr owave 
Hinshaw etHl 12008: 



Gold et al.1 120081: iNolta et all 12008 : iDunklev et al.1 



2008; 



Komatsu et al . 2008) contain the most valuable cosmolog- 



ical informati on to study the lar ge -scale properties of t he 
universe fseelBennett et al 112003 A iHinshawet al.ll2003al lbL 



120071 : I Jarosik et al.1 120071 : ISpergel et al.1 l2007h for previous 
releases of WMAP data). One of these features concerns 
the hypothesis that the set of temperature fluctuations of 
the cosmic microwave background radiation (CMB) is a 
stochastic realization of a random field, meaning that its 
angular distribution on the celestial sphere is statistically 
isotropic at all angular scales. Examination of the three- year 
WMAP data confirms highly significant departures from 
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statistical isotropy at large angular scales (lAbramo et al.l 
20061 : lHutere3l2006l: Iwiaux et alJliooS IVielva et al.ll2 006. 



20071 : lLand fe Magueiiol 120071: IPark 1 P ark fe Gott II I 
2007 : ICopi et al.1 120071 : lEriksen et al.l 12007 
20071 ) , previously found also in first-year 



Bernui et al 



data 



WMAP 
JTegmark. de Oliveira-Costa fe Hamilton! 120031: 
de Oliveira-Costa et al.l |2004| ; lHansen, Bandav &; Go rski 
120041: 



Hansen et al. 



2004 

2004: Biclcwicz 

2004al. 



2005 



ICopi 
Gorski fc Bandayl 



Land & Magueiio 



Huterer &; Starkman 

Erikse n et al. 

ICopi et al 




120061 ; iBernui et al.l I2006T ). Evidences for such large- angle 
anisotropy come from the asymmetry of the CMB angular 
correlations between the northern and southern ecliptic 
hemispheres (hereafter NS-asymmetry), with indications 
of a pr eferred axis of maximum hemis p herical asym- 
metr y (IHansen. Bandav &; Gorskil 
iBielewicz, Gorski &; Banday 




20051: lLand fc Magueiiol 12005a 




Hanse n et al 
Erikse n et al 
IBernui et al 



20061 : lEriksen et al.1 120071 : IBernui et al.1 120071 ). Further 
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more, other manifestations of large-angle anisotropy 
include the unlikely quadrupole-octopole planes alignment 
(refering both to the strong planarity of these multi- 
poles as w ell as to the alignment between such p lanes , 
see, e.g., iTegmark, de Oliveira-Costa fc Hamilton 12003: 



de Oliveira-Costa et al.1 12 004: Weeks 2004; 



Copi et al 



2006; IWiaux et all 120061: lAbramo et al.1 120061 : ICopi et al 
20071 : [Park. Park fe Gott IlJ I2007T ). Actually, sev- 
eral anomalous results were alre ady reported after 



analyses of the W MAP data (jChiang et al.l 



Vielva et al.l 



2004; ICopi, Huterer & Starkmanl 



2003 



2004; 



Bielewicz et al.ll2005l: ICruz et al.|l2005l. 120071: iMcEwen et all 
20061 : iBernui. Tsallis fe Villelal l200d. I2007T). For 



differ ent poi nt of view see, e.g., lHaiian fc Souradeepl 
(I2003T) 



Haiian, Souradeep & Cornish 



, ....... ,200^ 

Donoghue & Donoghud (120051 ) ; I Souradeep fe Haiianl (p004j 
Haiian & Souradeep (2005); Souradeep & Haiian (2005) 
Sour adeep, Haiian &; Basakl ([20061 ); lHaiian &; Souradeepl 



((200 

Possible sources of these 
elude non-C MB conta m inants, 
groun d s (lEriksen et al.l l2004bl: 



2005; de Oliveira-Costa fc Tegmarkl 



anomalies 

l ike residual 

IWibii 



m- 
fore- 



2006 



& Wolfendale 
ICruz et al 



2006; Abramo, Sodre & WuenscheT 120061 : IChiang et al 



2007; Lopez-Corredoira 2007) 



incorrectly subtracted 
dipol e and /or dynamic quadrupole ter ms ()Copi et al.l 



120061 : iHelling et all 120061: ICopi et al.1 120071), or systematic 
errors ([Schwarz et all 120041 : brain et al.1 l2007h . For this, 
particular efforts have been done by the WMAP team 
in last releases in order to improve the data process- 
ing by minimizing the effects of foregrounds (mainly 
coming from diffuse Galactic emission and astrophysical 
point-sources), artifacts (in the mapmaking process, in 
the in strument c haracterizati on, etc.), and systematic 
error s (jjarosik et al.ll2007l : Iffinshaw et alll2008l : iGold et all 
2008). As a result, the data released by the WMAP team 
include the Internal Linear Combination (hereafter ILC- 
5yr) full-sky CMB m ap suitable for large-angle temperature 
fluct uations studies (|Hinshaw et al.ll2007l . 120081 : iNolta et all 
2008). Here we investigate the ILC-5yr map, and for 
completeness, we also study the other full-sky cleaned 
CMB maps that were differently processed from WMAP 
five- and three-year data releases in order to account for 
foregrounds and systematics . Thus, we also consider the 
Kim- Naselskv-Christense n (|Kim, Naselskv &; Christensenl 
l2008h . W MAP-3vr ILC dHinshaw et al.1 20071) de Oliveira- 
Tegmark (Ide Oliveira-Costa &; Tegmark J2006I ) , and Park- 
Park- Gott |ParkIParkTGotniinEQ^ CMB maps, 
hereafter termed the HILC-5yr, ILC-3yr, OT-3yr, and 
PPG-3yr, respectively. 

A number of studies has been done looking for 
physical explanations of the above mentioned anomalies, 
specially searching for a unifying mechanism relating 
both the NS-asymm etry and the CMB lower multipoles 
align m ent (see, e.g., Wiaux et al. 2006; iRakic &; Schwarj 
120071 : iDvorkin. Peiris fe Hul 120071 ). Additionally, some 
processes that breaks down statistical isotro py during the 
inflat ionary epoch have been suggested (| Gordon et al.l 
20051: iGumrukcuoglu. Contaldi fe Pelosol 120071: 



Pulle n &; Kamionkowsil 120071 : lAckerman, Carrol &; Wisel 
20071 ) . Nonetheless, one can also interpret such CMB 



in this sense globally axisymmetric space-times have 
been proposed to accou nt for the mapped pr e ferred - 



axis in WMAP maps (Aurich 1 Lustin &; Steinerl 2005 

" 20051 : 



Hipolito-Ricaldi & Gomero 



2006: 



Land & Magueij 



Cresswell etHI l200d : 



_ Jaffe et al.1 120061: 
Gosh, Haiian &; Souradeepl 120071 ). Previous works in- 
vestigated the effects of a primordial homogeneous 
magnetic field on the CMB temperature fluctu ations on 
all a n gular scales (see, e.g., | Durrer, Kahniashvili fc Yatesl 



1 19981 : iDemiahski & Doroshkevichl l2007h . Here we study 
such primordial field as a possible physical mechanism 
to produce two large-angle phenomena, that is, the CMB 
NS-asymmetry and the lower CMB multipoles alignment. 

In section [2] we present this primordial magnetic field 
scenario, and show the effect of such a field on the CMB 
temperature fluctuations. Then, in order to reveal such 
effects on simulated CMB maps we develop a geometrical- 
statistical method, which is presented in section [3] After 
that we use our anisotropic indicator to perform, in sec- 
tion 21 the analyses of both sets of data, the Monte Carlo 
simulated CMB maps as well as the WMAP maps. At the 
end, in section \5\ we discuss our results and formulate our 
conclusions. 



2 PRIMORDIAL MAGNETIC FIELDS 
SCENARIO 

There is strong observational evidence for the pres- 
ence of large-scale intergalactic magnetic fields of fews 
/xG, and magnetic fields of similar strength within 
clust ers of galaxi es (see, e.g., iKrause, Beck &; Hum mel 



1989; 



Wolfe, Lanze tta &: Ore: 



n 



Clarke, Kroen berg &; Boehringerl l200ll : IWidrowl 



1992 



2002 



Xu et al.1 120061 ) . Nowadays it is believed that these mag- 



netic fields are amplifications of small primordial magnetic 
fields of the order of few nanoGauss, that would have 
occured due to different proc e sses, like galactical dy- 
namo (see, e.g., IParkerl Il971: |Vainshtein &; Ruzma ikin 
Il972l : IVainshtein &; Zel'dovichl Il972l ), during anisotropic 
protogalactic colla pses, or due to d i fferential rotation in 
galaxies (see, e.g., Piddin gtorj Il970l : iKulsrud &; Anderson! 
1992). In turn, such primordial seeds of magnetic fields 
would have seve ral origins, for instance an electroweak 
phase transition (| Vachaspati| ll993l : iKibble &; Vilekinl Il995l : 
iBavm, Bodeker fc McLerranl 19961) or quark-hadron phase 



large-angle anisotropy as being of cosmological origin, 



transition (see, e.g.. lOuashnock, Loeb fc Spergel [l989; 
Che ng &; Olintol Il994). Here, we assume the existence of 
a primordial homogeneus magnetic field and investigate 
their effects on the CMB temperature fluctuations at 
large- angles. 

According to the Einstein equations for linearized met- 
ric perturbations, in absence of a magnetic field, the vector 
metric perturbations go like a~ 2 and the velocity induced by 
these perturbations goes like a -1 (where a is the scale fac- 
tor that accounts for the expansion of the universe) . There- 
fore, t hey decay very fast with the expansion of the uni- 
verse (|Mukhanovl 12005) and the velocities produced by vec- 
tor metric perturbations do not contribute significatively to 
the CMB temperature fluctuations. 

The presence of a homogeneus magnetic field in the 
early universe change this scenario because such fields mod- 
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ify the behavior of charged particles in the primordial plasma 
via Lorentz forces producing additional velocity gradients 
in the fluid. In this way, a magnetic field induces Alfven 
waves in the primordial plasma that propagate at velocity 
vj± , changing the speed of soun d in the photon-bary on fluid where 



;on-bary c 

2 c 2 s + v\ cos 2 0, where f Adams et al.lll996h 



2 

V A 



47r(p + p) 



(1) 



with p and p being the density and pressure in the ra- 
diation dominated era, respectively, Bo is the strength of 
the magnetic field B, is the angle between B and the k- 
mode of the Fourier expansion of Alfven velocity vj±. These 
Alfven wave modes induce small rotational velocity pertur- 
bations which, for the scales of our interest, have the form 
v ~ vovj^kt cos #, where t is the cosmic time, k = \k\, and 
vo is the initial velocit y, which we assume to have a power 
spectrum of the form (|Durrer, Kahniashvili &; Yatesl Il998l ; 
IChen et al.|[20oi ) 



(voi vqj) oc k n (8ij - h kj) . 



(2) 



Vectorial contributions to CMB temperature 
fluctuati ons are present via Doppler a nd Sachs- Wolfe 
effects ([Durrer, Kahniashvili fc Yateslll998T ) 



f (n) (vec) = 



-Vt • n 



to 



f 



V • hdt , 



(3) 



where h indicates any direction in the sky, Vt = V -\- 
v, v is the vectorial metric perturbation, V is the ve- 
locity produced by z7, and the integration is between 
the actual time to and the decoupling time tdec- Thus, 
only the rotational velocity perturbations v contributes 
effectively to temperature fluctuations because v and V de- 
cay q uickly with the expansion of the universe (|Mukhanovl 
120051 ). Therefore, the signatures of a homogeneus primor- 
dial magnetic field B on the CMB temp erature fluctuations 
are (JDurrer, Kahniashvili fc Yateslll998h 



6T 1 



(n, k) ~ h • vo va k tdec cos . 



(4) 



As we can see, the existence of a preferred direction as- 
sociated to B affects the CMB temperature fluctuations 
through the angle #, which implies, a straight dependence 
on the orientation of the vector B, consequently, a break 
down of the statistical isotropy in the CMB sky. 

The most suitable quantity to simulate statistically 
anisotropic CMB skies is the correlation matrix of multi- 
pole coefficients. For this, we expand the sky temperature 
fluctuations in spherical harmonics 

ST, 



-(h) 



Cl£ m Y£m(h) . 



a£ri 



5T 



(n)Y£ m (h)dQ, 



(5) 



(6) 



and then use eqs. (|2]), Q and (0 to calculate the correlation 
matrix {a£ rn a\, rn ,) . In this work we assume a scale- invariant 
power spectrum for the magnetic field, i.e., n — —5 in the 
eq. ©. 

The explicit calculation of the correlation 
matrix elements, for a Harrison- Zel'dovich scale- 
invariant power spectrum for the magnetic field, 



gives (IDurrer, K ahniashvi li &; Yatesl Il998l ; IChen et al.l 
I2OO4I ; iNaselskv fe Kimll2008h 

{aimCLl'm') = &mm'[&U'Ci rn + ((^ + i^/_i + ^_i^/ + i)^ m ] , (7) 



Cf m = 2,12 x 10- (*, ) x 

2l 4 + 4£ 3 - f - 31 + 6m 2 - 2£m 2 - 2£ 2 m 2 

(£ + 2)(£ + 1)£(£ - 1)(2^ - 1)(2^ + 3) 
and 



(8) 



U £m 


9tt ^(£ + m + l)(£-m + l){£ + m){£-m) 




32 2£ 4 + 4P - P - 3£ + 6m 2 - 2£m 2 - 2£ 2 m 2 




y/(2£-l)(2£ + S)(£ -!)(£ + 2) 



(2£ + 1) 

noticing that Bo is given in nanoGauss (nG). Thus, the com- 
putation of the correlation matrix leads to non-zero elements 
of the type (£, m) = (£! , vn') and (£, m) = {£' ± 2, m), while 
all the other elements are zero. In other words, the multipole 
matrix correlation is non-diagonal which is an inheritance of 
its angular dependence on and vo, shown in the eq. (HJ), 
and this fact being a consequence of the presence of the 
magnetic field at early times. In this scenario, correlations 
appear between different scales, like £ and £ ± 2, and then 
the multipole coefficients a^ m at such scales are correlated. 
For this reason, one concludes that large-angle anisotropy 
features should appear in the CMB maps produced by this 
primordial magnetic field. If we consider a magnetic field in 
a ACDM universe the correlation matrix will be 



(9) 



where 



{a£mCl£>', 



\ACDM 



% d£l. 



B 



* \ACDM _ ^ACDM x x 

£md£'m / ) ~ ^ '£ ^mm' °l 



(10) 



(ii) 



We observe that correlations between temperature 
fluctuations do not only depend on the angular sepa- 
ration between two points, but also on their orientation 
with respect to the magnetic field, that is, m- dependence. 
This m-dependence, caused by the fact that Bo / 0, means 
that now the CMB temperature fluctuations are actually 
statistically anisotropic. 

Our purpose here is just to illustrate the effect in- 
duced by these large-angle anisotropics in CMB maps, 
and compare them with those features found in the CMB 
maps from WMAP. To study the relationship between 
low-^ (i.e., from £ — 2 to £ — 10) anomalies in such 
anisotropic magnetic field scenario, we produce, according 
to eqs. (ff))- (pT]) . five sets (with different strengths Bo) of 
Monte Carlo simulated CMB maps. For the statistically 
isotropic part (that is, the first term in eq. (|10J)) we use the 
CMBFAST tool to calculate the angular power spectrum 
C^ , and after that we use Cholesky decomposition of 
the matrix (|10[) in order to deal with the non-diagonality. 
Then we randomically simulate the a£ m coefficient sets 
and from them we generate the CMB temperature 
fluctuations maps. The simulations were performed for 
magnetic fields intensities Bo = (20 ± 10) nG, to be in 
agreement with known limits at co smological scales for the 
n = —5 power spect ru m (see , e. g., iBarrow, Ferreira &: Silk 
Il997l : JChen et al.l 12004 INaselskv fe Kiml 120081 . 
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and iKahniashvili. Maravin &; Kosowskvl 120081 for other 
power spectrum indices). 



statistical isotropy hypothesis. As we shall see, the sigma- 
map analysis is able to reveal large- angle anisotropics such 
as the NS-asymmetry. 



3 THE 2PACF AND THE SIGMA-MAP 
METHOD 

Our method to investigate the large-scale angular correla- 
tions in CMB temperature fluctuations maps consists in 
the compu tation of the 2-point angular correlation function 
(2PACF) ([Padmanabhan 1993) in a set of spherical caps 
covering the celestial sphere. 

Let Q 70 = Q(6jj(j)j] 70) C <S 2 be a spherical cap region 
on the celestial sphere, of 70 degrees of aperture, with ver- 
tex at the J-th pixel, J — 1, . . . , A^caps, where (6j, <j>j) are 
the angular coordinates of the J-th pixel's center. Both, the 
number of spherical caps A^caps and the coordinates of their 
centers { Oj,<j)j) are defined using the HEALPix pixelization 
scheme (|G6rski et al.ll2003r ). The union of the Acaps spher- 
ical caps covers completely the celestial sphere S 2 . 

Given a pixelized CMB map, the 2PACF of the temper- 
ature fluctuations ST corresponding to the pixels located in 
the spherical cap Qj is defined by (|Padmanabhanl 1 1993T ) 



C( 7 r ={5T{^^)5T{0i',4>i'))- 



(12) 



where cos 7 = cos Oi cos 6^ +sin Oi sin 0^ cos(0i— fat ) , and 7 £ 
(0, 270] is the angular distance between the i-th and the i'-th 
pixels centers. The average ( ) in the above equation is done 
over all the products 6T(9i,fa)6T(0i/,fat) such that 7^ = 
7 E (0 - 1)6, kS], for k = 1, ..., AT bins , where 6 = 27o/AT bins 
is the bin- width. We denote by C^ = C(7fc) J the value of 
the 2PACF for the angular distances 7^ G ((k — 1)6, kS]. 
Define now the scalar function a : Q^ C <S 2 1— ► 5R + , for 
J = 1, . . . , A^caps, which assigns to the J-cap, centered at 
(6j,(j)j), a real positive number aj = cr(0j,(j)j) 6 5R + . The 
most natural way of defini ng a measure a is through the 
variance of the C^ function (jBernui et al.ll2007l h 



bins 



2 



ATt 



bins 



£ ( C * J ) 



(13) 



To obtain a quantitative measure of the angular correla- 
tions in a CMB map, we cover the celestial sphere with 
Acaps spherical caps, and calculate the set of sigma values 
{aj, J = 1, ..., Acaps} using the eq. ([T3]) . Associating the 
J-th sigma value aj to the J-th pixel, for J = 1, ... , Acaps, 
one fills the celestial sphere with positive real numbers, 
and according to a linear scale (where a miminiain _> blue, 
^maximum _^ rec ^ one conver ts this numbered map into 
a coloured map: this is the sigma-map. Finally, we find 
the multipole components of a sigma-map by calculat- 
ing its angular power spectrum. In fact, given a sigma- 
map one can expand a — a(0,<f>) in spherical harmon- 
ics: (t{0,(J)) — ^2 £ Ai mYi m(0 , fa) . Then the set of values 

{S t ,l = 1,2,...}, where S e = (l/(2£ + 1)) Y? m =-z \ A ^\^ 
give the angular power spectrum of the sigma-map. 

A power spectrum S\ of a sigma-map computed 

from a given WMAP map, provides quantitative information 
about large-angle anisotropy features of such a CMB map 
when compared with the mean of sigma-map power spectra 
obtained from Monte Carlo CMB maps produced under the 



4 DATA ANALYSES AND RESULTS 

Now we shall apply the sigma-map method to scrutiny the 
large- angle correlations present in the WMAP maps, and 
perform a quantitative comparison with the result of a sim- 
ilar analysis performed in sets of simulated Monte Carlo 
(MC) CMB maps. These sets of simulated maps are pro- 
duced according to the primordial magnetic field scenario 
discussed above, where we consider five cases, for an equal 
number of values for the parameter Bo, to examine the 
influence of the field intensity on the CMB angular correla- 
tions. These analyses let us to test the hypothesis that the 
anomalous angular correlations found in the WMAP data 
could be explained by the presence of a magnetic field act- 
ing in the decoupling era. The WMAP data under inves- 
tigation are the full-sky cleaned CMB maps derived from 
WMAP five- and three-year releases, namely the ILC-5yr, 
HILC-5yr, ILC-3yr, OT-3yr, and PPG-3yr CMB maps. Be- 
cause we are interested in understand the possible cause- 
effect relationship between a primordial magnetic field and 
the large- angle anomalies mapped in CMB WMAP data, 
we concentrate our study on the low-order multipoles range 
£ = 2-10. 

For our analyses we produce five sets of 1 000 MC 
CMB maps each, corresponding to the cases when Bo = 
0, 10, 15, 20, and 30 nG. The case Bo = means that the 
MC maps were g enerated using a pure ACDM angular powe r 
spectrum seed (|Spergel et al.1 120071 : iKomatsu et al.1 l2008h . 
and this case refers to the statistically isotropic CMB maps. 

For Bo y^ the simulated MC CMB maps were pro- 
duced considering the two contributions to the random 
a£ m -modes according to eq. ([TO]) , that is, the statistically 
isotropic part plus the component due to the magnetic field. 
As a result of this, the mean angular power spectra of each 
set of MC maps, for Bo = 10, 15, 20, 30 nG, satisfies the 
Sachs- Wolfe plateau effect. However, due to the contribu- 
tions (aimCLgrmf) (see the second term in eq. (|10[0 the 
plateau of these mean angular power spectra are shifted up- 
wards, where the value of such shifts is proportional to Bo. 
In order to compare features obtained, through the sigma- 
map method, from MC maps with power spectra that be 
consistent with the ACDM power spectrum we normalize 
the MC power spectra corresponding to the Bo ^ cases. 
The normalized mean angular power spectra Gg of these four 
sets of MC maps Bo / cases and the angular power spec- 
trum of the ACDM model, plus its cosmic variance limits, 
are shown in figure [1] 

In our simulations we have assumed that the magnetic 
field is pointing in the South Galactic Pole-North Galactic 
Pole (SGP-NGP) direction. Notice that all the sky maps 
plotted here are in galactic coordinates, which means that 
the equator of the map corresponds to the Galactic plane, 
and the axis SGP-NGP is perpendicular to this plane. 

The data analyses consist on the following steps. For 
each MC temperature map we compute its corresponding 
sigma-map (hereafter called sigma-map-MC), and then we 
calculate its angular power spectrum {Si, I = 1,2, ...}. The 
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Figure 1. Normalized mean angular power spectra of the MC 
CMB maps used to produce the sigma-maps-MC. The bullet, 
plus, triangle, square, and X symbols represent the data when 
Bo = 0, 10, 15, 20, 30 nG, respectively. The dashed lines are the 
cosmic variance limits of the ACDM (i.e., Bo = 0) power spec- 
trum case. Note that the four sets of MC maps with Bo 7^ have 
their angular power spectra consistent with the ACDM case. 



statistical significance of the sigma-map angular power spec- 
tra computed from WMAP data (hereafter called sigma- 
map- WMAP) comes from the comparison with the five sets 
of 1 000 angular power spectra obtained from the sigma- 
maps-MC. To illustrate the effect of the magnetic field in 
the CMB low-order multipoles in MC maps, we show three 
cases in figure [2] the mean of 100 sigma-maps-MC obtained 
from a similar number of MC computed considering the sta- 
tistically isotropic case Bo = (top), and considering mag- 
netic field intensities Bo = 10 nG (middle) and Bo = 30 nG 
(bottom) respectively. In figure [3] instead, we show three 
sigma-maps with large value of the dipole term Si , indica- 
tive of the hemispherical asymmetry phenomenon, produced 
from the ILC-5yr map (top) and from MC CMB maps with 
different magnetic field intensities: Bo = 10 nG (middle) 
and Bo = 30 nG (bottom). 

In figure S] we present the results of the sigma-maps 
spectra analyses. We observe that the sigma-maps- WMAP, 
obtained from the ILC-5yr, HILC-5yr, ILC-3yr, OT-3yr, and 
PPG-3yr CMB maps, reveal a dipole moment S 1 larger 



than 95% of the values S ± 



MC-ACDM 



corresponding to sigma- 



maps-MC obtained from statistically isotropic CMB maps. 
This fact indicates that the NS-asymmetry phenomenon is 
present in WMAP data at 95% CL. However, we also ob- 
serve in figure [4]t hat the angular power spectra of the sigma- 
maps-MC varies with the magnetic field intensity, thus to 
larger values of Bo correspond sigma-maps-MC with larger 
values of the terms Si (dipole), & (quadrupole) , S3 (oc- 
topole), etc. In other words, a larger value of Bo produces a 
stronger hemispherical asymmetry in the MC CMB maps. 
Therefore, we conclude that for sufficiently large Bo, one 
can interpret the spectra resulting from sigma-maps- WMAP 
as being not anomalous at all, but consistent with those 
CMB temperature maps produced according to the primor- 
dial magnetic field scenario. Interestingly, simulations seems 




Figure 2. These sky maps, in Galactic coordinates and from top 
to bottom, represent the mean of 100 sigma-maps-MC obtained 
from a similar number of MC computed considering three cases: 
Bo = 0, Bo = 10 nG, and Bo = 30 nG, respectively, and with 
the magnetic field pointing in the SGP-NGP direction, which 
means that the equator is the preferred plane. Notice that for 
Bo / the region around the equator concentrates strong angular 
correlations, here represented by the intense and large red sky 
patches. 



to indicate that an intensity of Bo ~ 15 nG is enough to have 
a NS-asymmetry as strong as WMAP data (see figure H}. 

The effect of the magnetic field on the set of (£, m)- 
modes, for a given multipole £, deserves a close inspection. 
We know that in the statistically isotropic case the power 
of the ^-multipole in a CMB map, given by Ct — (l/2£ + 
1) Ylm=-e \ a ^ m \ 2 5 * s uniformly distributed between the (2£-\- 
1) modes, for m = —£, ...0, ...£. However, when a primordial 
magnetic field is acting on the CMB, and pointing in the 
SGP-NGP direction, we expect a non-uniform distribution 
of such power with a net predilection for the planar modes 
perpendicular to this axis, event that is termed the m = £- 
preference. 

To investigate if the multipoles of our MC are experi- 
encing this planarity phenomenon when Bo 7^ 0, we perform 
a probabilistic analysis of how the power of the quadrupole 
moment C2 and octopole moment C3 are distributed in their 
(21 + 1) = 5 and (2£ + 1) = 7 modes, respectively. For this, 
as a criterium for measuring the predilection for the (£, ni)- 
mode we consider those values obtained for the quadrupole 
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Figure 3. For illustration we show three sigma-maps: the sigma- 
map-WMAP from the ILC-5yr CMB map (top), and two sigma- 
maps-MC having large dipole term Si one obtained from a MC 
map with Bo = 10 nG (middle) and the other obtained from a 
MC map with B = 30 nG (bottom). 



and octopole from the ILC-5yr CMB map, that is, (£, m) is 
a preferred mode when it takes more than 35% of the power 
of such multipole, i.e., \a im \ 2 + \a} m \ 2 > 0.35(2^ + \)Cg. 
Given a subset of MC where the power of at least one of 
their (2£' + l)-modes, of a given £'- multipole, is greater 
than 0.35 (2£ f + 1)CV, we define Vi'm' as the probability 
that the mode (£ / ,m / ) satisfies the power distribution cri- 



terium (PDC): \ay. 



+ \or. 



> 0.35 (2£' + \)C V . For 



instance, ^2,2 is the probability that the quadrupolar mode 
(£',m f ) = (2,2) satisfies |a 2 , 2 | 2 + \a^ 2 \ 2 > (0.35) (5) C 2 , 
where such probability is computed considering the subset of 
MC where at least one of the 5 modes (£\ ra'), for £' — 2 and 
m — —2,-1,0,1,2, has power greater than (0.35) (5) C2. 
We perform a comparative analysis for two sets of MC data, 
namely the set of 1000 MC ACDM (i.e., B = 0, here- 
after MC-ACDM) and the set of 1 000 MC produced with 
Bo = 20 nG (hereafter MC-B20). 

Additionally, to realize a possible correlation between 
the quadrupole-octopole planes alignment and the NS- 
asymmetry phenomena we also investigate the occurrence 
of the m — ^-preference in three subsets of the MC-ACDM 
and MC-B20 sets, namely those subsets that satisfies the 
PDC for the quadrupole and the octopole simultaneously 
and such that their corresponding sigma-maps-MC have: (i) 




3 4 5 
multipoles (I 



Figure 4. Angular power spectra of the sigma-maps-MC for MC 
produced with different magnetic fields intensities together with 
the sigma-maps- WMAP spectra. The plus, square, X , bullet, and 
triangle symbols represent the data from the ILC-5yr, HILC-5yr, 
ILC-3yr, PPG-3yr, and OT-3yr CMB maps, respectively. From 
bottom to top, the solid (dashed) lines represent the mean (95% 
CL) values computed from the set of sigma-maps-MC correspond- 
ing to the cases Bq = 0, 10, 15, 20, 30 nG, respectively. 



any dipole moment value Si , hereafter denoted by MC- 
ACDM 0-SD and MC-B20 0-SD subsets, respectively; (ii) 
dipole moment value Si larger than the mean value plus 
one standard deviation, hereafter these data are termed MC- 
ACDM 1-SD and MC-B20 1-SD, respectively; (hi) dipole 
moment value Si larger than the mean value plus two 
standard deviations, hereafter these data are termed MC- 
ACDM 2-SD and MC-B20 2-SD, respectively. The fraction 
of MC, in the MC-ACDM and MC-B20 sets, that satisfies 
the above PDC for the quadrupole and the octopole simul- 
taneously is ~ 0.63 for all these data subsets. 

Our results for these computations are shown in Table [1] 
where we observe the following results. First, as expected for 
the three sets of MC-ACDM data, the analysis reveals that 
the quadrupole and the octopole have their power uniformly 
distributed between all the (£, ra)-modes. Second, regarding 
the MC-B20 data sets, it is observed a weak preference for 
the (£, m) 7^ (-£, 0) modes in 0-SD and 1-SD data sets, while 
the planarity or m = ^-preference is actually evident in the 
set 2-SD. Due to this fact, and according to the definition 
of Vim, we conclude that there is a net correlation between 
the quadrupole-octopole planes alignment (represented by 
the m — ^-preference through the large values of 7^2,2 and 
7^3,3) and the NS-asymmetry phenomena (represented by 
the fact that these large probability values appear only con- 
sidering those MC that produce sigma-maps-MC with the 



largest dipoles £1 ) . We understand this m — ^-preference 
as a consequence of the planarity induced by the preferred 
direction settled by the magnetic field, as illustrated with 
two examples in figure [5] Additionally, one notices the effect 
that this planarity produces in the sigma-maps-MC causing 
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Figure 5. The left (right) panels are the quadrupole £ = 2, the 
octopole £ = 3 and their sum £ = 2 + 3 corresponding to the 
MC CMB map that produces the sigma-map-MC plotted in the 
middle (bottom) panel of figure [3] 



a concentration of strong angular correlations (red regions) 
around the equator, which is the preferred plane, as clearly 
seen in the mean of the sigma-maps-MC skies showed in 
figure [2] 

A important part of our analyses are the robustness 
tests. To assert the robustness of our results when using a 
different set of parameters in the sigma-map method, we 
investigated the effect of changing 70, Afk ms , iV C aps, and 
^side °f t ne CMB map in analysis, in the computation of the 
sigma-maps. For this we performed several sigma-map cal- 
culations using spherical caps of 70 = 45°, 60° of aperture, 
^bins = 45 > 60 > 90 > and ^caps = 768, 3 072, resulting in mi- 
nor differences in all these cases. Additionally, we also exam- 
ine the influence of the angular resolution of the CMB maps 
by computing the sigma-maps considering different pixeli- 
sation parameters of the CMB maps, namely N s - l( ^ e = 16 
and 32. In particular, the sigma-maps showed in figure [2] 
plotted in galactic coordinates, and their corresponding an- 
gular power spectra analyses in figure [1] were calculated 
using 70 = 45°, A^bins = 45 , ^Vcaps = 768, and iV side = 32. 
Summarizing, our robustness tests show that using a set of 
parameters within a certain range of values, we obtain re- 
sults that are fully consistent with those showed in figure 2] 



5 CONCLUSIONS 

We investigated a plausible primordial scenario where a ho- 
mogeneous magnetic field acting on the photon-baryon fluid 
at the recombination era, introduces a preferred direction 



that establish a statistical isotropy breaking. This setting 
offers a possible explanation for those anomalies found in 
WMAP data that seems to be associated to a preferred axis 
in the space, like the hemispherical NS-asymmetry and the 
planarity of some CMB low-order multipoles (where such a 
plane is perpendicular to that axis). It is found that this 
magnetic field induces correlations between the CMB mul- 
tipoles, manifested through non-diagonal terms in the mul- 
tipole correlation matrix (see eqs. Q-JSJ)). Accordingly, we 
have simulated five sets of CMB skies considering an equal 
number of strengths Bo of the magnetic field. These MC 
CMB maps, with multipoles in the range £ = 2 — 10, are used 
to investigate their large-angle correlations using our sigma- 
map method. With these data sets we performed a quanti- 
tative analysis of their large-angle signatures and perform 
a comparison with similar features computed from WMAP 
maps. For this, we consider a set of full-sky cleaned CMB 
maps, obtained from five-year and three- year WMAP data 
releases, namely the ILC-5yr, KCN-5yr, ILC-3yr, OT-3yr, 
and PPG-3yr maps, which were differently processed in or- 
der to account for foregrounds and systematic errors. 

Our results can be summarized as follows. First, our 
analyses corroborate that an uneven hemispherical distri- 
bution in the power of the large-angular correlations, best 
known as NS-asymmetry phenomenon, is present in all these 
WMAP maps at more than 95% CL as compared with sta- 
tistically isotropic CMB maps produced according to the 
ACDM cosmological model (i.e., Bo =0). Second, we found 
that this hemispherical asymmetry phenomenon, is such 
that higher is the value of the field strenght .Bo, greater 
is the probability that such asymmetry be present in the 
MC CMB maps, and be revealed with a high dipole value 
S{ in its corresponding sigma-map-MC. In other words, 
our results show that the correlations introduced by a mag- 
netic field mechanism in low-order CMB multipoles, for 
a sufficiently intense field, reproduce the NS-asymmetry 
present in WMAP data as a common phenomena and not 
as an anomalous one. Third, as expected, the MC CMB 
maps with Bo 7^ exhibit the planarity effect in their low- 
order multipoles (see, for illustration, figure [5} as quan- 
tified in the Table [1] where the ^-multipole power is con- 
centrated in the m = ±£ modes, and as evidenced by the 
intense red spots representative of strong angular correla- 
tions appearing in the equatorial region of the sigma-maps- 
MC (see, for illustration, figure [2}. Fourth, as shown in the 
Table [T] we found a significant correlation between the NS- 
asymmetry and quadrupole-octopole (planarity and) align- 
ment phenomena. For completeness, we have also verified 
that our results are robust under different sets of parame- 
ters, as mentioned in the previous section, involved in the 
sigma-map method calculations. 

The large angular scales CMB anomalies challenge the 
statistical isotropy expected in the ACDM concordance 
model. Our results suggests that perhaps some of them could 
be suitable explained by just one physical phenomenon. As 
a matter of fact, there are several attempts to find the origin 
of large- angle CMB anomalies, but further investigations are 
needed to fully comprehend if they have one or more causes. 
In this sense, statistically anisotropic models explaining sev- 
eral anomalies with a minimum set of hipothesis should be 
explored. 

We believe that our results shall motivate the study 
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Table 1. Probability comparative analyses of the power distribution for the a£ m -modes of the quadrupole (£ = 2) and the octopole (i = 3) 
in those subsets of MC-ACDM and MC-B20 that satisfies the PDC for the quadrupole and octopole simultaneously. The probability 
Vim is defined in the text. 0-SD means that we are consider for analysis the subset of the MC sets (mentioned above) such that their 
corresponding sigma-maps-MC have any dipole moment value S^ c ; 1-SD (2-SD) means that we are consider for analysis the subset of 
the MC sets (mentioned above) such that their corresponding sigma-maps-MC have dipole moment value S^ C larger than the mean 
value plus one (two) standard deviation(s). 



MC data \ Ve„ 



V 2 ,o 



T>2, 



V 2 ,: 



V 3 ,o 



T>3, 



V 3 ' 



V 3 ,: 



MC-ACDM 0-SD 


26.7% 


39.0% 


34.3% 


26.3% 


23.8% 


26.6% 


23.3% 


MC-ACDM 1-SD 


25.8% 


39.9% 


34.3% 


25.9% 


22.4% 


27.9% 


23.8% 


MC-ACDM 2-SD 


30.2% 


32.6% 


37.2% 


20.2% 


26.6% 


26.6% 


26.6% 


MC-B20 0-SD 


10.4% 


45.0% 


44.6% 


13.0% 


29.7% 


27.6% 


29.7% 


MC-B20 1-SD 


11.4% 


39.3% 


49.3% 


7.5% 


31.9% 


25.3% 


35.3% 


MC-B20 2-SD 


4.0% 


32.0% 


64.0% 


0.0% 


14.3% 


28.6% 


57.1% 



of other statistically anisotropic scenarios, in particular, 
those where correlations between CMB anomalies ap- 
pear. The magnetic field scenario analysed here is just 
the simplest one, and considering more realistic fields 
new effects could be present on the CMB tem perature 
and polarization da t a (see , e.g., |G iovannin il 120061 2007; 
Giovannini fe Kunzd 120071: iKahniashvili fe Ratral 12 005 



Subramaniam Seshadri fc Barrow 1120031; iKosowskv &; Loebl 
19961 ; ISeshadri &; Subramanianl l200lL and references 
therein). In the same form, as well, possible relationships 
between these new effects and statistically anisotropic 
phenomena found in WMAP data could be established, all 
this deserving a better investigation. 
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